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In this paper we study state-space realizations of Linear and Time-Invariant (LTI) systems. Mo- 
tivated by biochemical reaction networks, Goncalves and Warnick have recently introduced the notion 
of a Dynamical Structure Functions (DSF), a particular factorization of the system's transfer function 
matrix that elucidates the interconnection structure in dependencies between manifest variables. We 
build onto this work by showing an intrinsic connection between a DSF and certain sparse left coprime 



factorizations. By establishing this link, we provide an interesting systems theoretic interpretation of 
sparsity patterns of coprime factors. In particular we show how the sparsity of these coprime factors 



allows for a given LTI system to be implemented as a network of LTI sub-systems. We examine possible 



applications in distributed control such as the design of a LTI controller that can be implemented over 
a network with a pre-specified topology. 
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I. Introduction 

Distributed and decentralized control of LTI systems has been a topic of intense research 
focus in control theory for more than 40 years. Pioneering work includes includes that of Radner 
JD, who revealed the sufficient conditions under which the minimal quadratic cost for a linear 
system can be achieved by a linear controller. Ho and Chu [0, laid the foundation of team theory 
by introducing a general class of distributed structures, dubbed partially nested, for which they 
showed the optimal LQG controller to be linear. More recently in [fTTTl . [[121 . [fT3l . Ifl4l important 
advances were made for the case where the decentralized nature of the problem is modeled as 
sparsity constraints on the input-output operator (the transfer function matrix) of the controller. 
These types of constraints are equivalent with computing the output feedback control law while 
having access to only partial measurements. Quite different from this scenario, in this work we 
are studying the meaning of sparsity constraints on the left coprime factors of the controller, 
which is not noticeable on its transfer function. In particular, we show how the sparsity of these 
coprime factors allows for the given LTI controller to be implemented over a LTI network with 
a pre-specified topology. 

More recently, network reconstruction of biochemical reaction networks have motivated a 
careful investigation into the nature of systems and the many interpretations of structure or 
sparsity structure one may define [18]. In this work, a novel partial structure representation 
for Linear Time-Invariant (LTI) systems, called the Dynamical Structure Function (DSF) was 
introduced. The DSF was shown to be a factorization of a system's transfer function that 
represented the open-loop causal dependencies among manifest variables, an interpretation of 
system structure dubbed the Signal Structure. 

A. An Introductory Example / Tiffl/ 

One important characteristic of the DSF is its ability to represent the impact that observed 
variables have on each other. This can often effectively describe the interconnection structure 
between component subsystems within a given system. Consider for example the 3-hop ring 
(also called "delta") network in Figured] where all the Q(s) and P(s) blocks represent transfer 
functions of continuous-time LTI systems. We denote with L(s) the transfer function from the 
input signals U(s) to the outputs Y(s). By directly inspecting the signal flow graph in Figure [T] 
we can write the algebraic equations: 
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Fig. 1. A 3-Hop Ring Network 
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We make the additional notation 

O O 
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and P(s) d = 
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(2) 



and we define ad-hoc the (Q(s), P(s)) pair to be the Dynamical Structure Function associated 
with the L(s) LTI system. (The rigorous definition of DFS will be introduced in Section HI 
following the original mathematical derivation from [ED.) An interesting observation, which is 
the main thesis of this work, is that the structure of the subsystems interconnections in Figure [T] 
is no longer recognizable from the input-output relation described by the transfer function of 
the aggregate system L(s) = (l 3 — Q{s)) X P{s) since the transfer function L(s) does not 
have any sparsity pattern and in general does not have any other particularities. The structure 
however, remains visible and it is captured in the quite particular sparsity patterns of Q(s) and 
P(s), respectively. This key property makes the DSF susceptible of becoming a perfectly suited 
theoretical concept to model any LTI network. 

We want to illustrate further how the DSF determines via equation (fT6l) the topology of the 
LTI network that can describe the given LTI system L(X). If we consider Qi3(s) identically zero 
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in (OQ) which would mean"breaking" the ring network from Figure Q] then it becomes a cascade 
connection and L(X) can be implemented as a "line" network. A "line" network controller could 
be interesting for example motion control of vehicles moving in a platoon formation. 

Note that, in general, the impact of observed variables on each other, represented by the DSF, 
and the interconnection between subsystems, can be quite different structures. This is because 
the states internal to one subsystem are always distinct from another, while the states internal to 
component systems in the DSF may be shared with other components. Nevertheless, the point 
in this example, that the DSF, as a factorization of a system's transfer function, captures an 
important notion of structure, is always true. Details about the distinctions between subsystem 
structure and the signal structure described by the DSF can be found in [[T9l . 

B. Motivation and Scope of Work 

In this paper we look at Dynamical Structure Functions from a control systems perspective. 
A long standing problem in control of LTI systems was synthesis of decentralized stabilizing 
controllers ([4]) which means imposing on the controller's transfer function matrix K(s) to have 
a diagonal sparsity pattern. Quite different to the decentralized paradigm, the ultimate goal of our 
research would be a systematic method of designing controllers that can be implemented as a LTI 
network with a pre-specified topology. This is equivalent with computing a stabilizing controller 
K(s) whose DSF (Q(s), P(s)) satisfies certain sparsity constraints [20] . So, instead of imposing 
sparsity constraints on the transfer function of the controller as it is the case in decentralized 
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control, we are interested in imposing the sparsity constraints on the controller's DSF. This 
would eventually lead to the possibility of designing controllers that can be implemented as a 
LTI network, see for example Figure |2] 

C. Contribution 

The contribution of this paper is the establishment of the intrinsic connections between the 
DSFs and the left coprime factorizations of a given transfer function and to give a systems 
theoretic meaning to sparsity patterns of coprime factors using DSFs. The importance of this is 
twofold. First, this is the most common scenario in control engineering practice (e.g. manufac- 
turing, chemical plants) that the given plant is made out of many interconnected sub-systems. 
The structure of this interconnection is captured by a DFS description of the plant which in turn 
might translate to left coprime factorization of the plant that features certain sparsity patterns on 
its factors. This sparisty might be used for the synthesis of a controller to be implemented over 
a LTI network. Conversely, in many applications it is desired that the stabilizing controller be 
implemented in a distributed manner, for instance as a LTI network with a pre-specified topology. 
This is equivalent to imposing certain sparsity constraints on the left coprime factorization of the 
controller (via the celebrated Youla parameterization). In order to fully exploit the power of the 
DSFs approach to tackle these types of problems, we find it useful to underline its links with the 
classical notions and results in control theory of LTI systems. We provide here a comprehensive 
exposition of the elemental connections between the Dynamical Structure Functions and the 
Coprime Factorizations of a given Linear Time-Invariant (LTI) system, thus opening the way 
between exploiting the structure of the plant via the DSF and employing the celebrated Youla 
parameterization for feedback output stabilization. 

D. Outline of the Paper 

In the second Section of the paper we give a brief outline of the theoretical concept of 
Dynamical Structure Functions as originally introduced in [Tl8l . In the third Section, we show 
that while the DSF representation of a given LTI system L(s) is in general never coprime, 
a closely related representation dubbed a viable (W, V) pair associated with L(s) is always 
coprime. We also provide the class of all viable (W, V) pairs associated with a given L(s). The 
fourth Section contains the main results of the paper and it makes a complete explanation of 
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the natural connections between the DSFs and the viable (W, V) pairs associated with a given 
L(s) and its left coprime factorizations. The last Section contains the conclusions and future 
research directions. In the Appendix A we have provided a short primer on realization theory 
for improper TFMs which is indispensable for the proofs of the main results. The proofs of the 
main results have been placed in Appendix B. 

II. Dynamical Structure Functions 

The main object of study here is a LTI system, which in the continuous-time case are described 
by the state equations 

x(t) = Ax(t) + Bu(t); x(t Q ) = x Q (3a) 
y(t) = Cx(t) + Du{t) (3b) 

where A, B , C \ D are n x n, n x m, p x n, p x m real matrices, respectively while n is also called 
the order of the realization. Given any n-dimensional state-space representation (l3at . (|3bl of 
a LTI system (A, B, C, D), its input-output representation is given by the Transfer Function 
Matrix (TFM) which is the p x m matrix with real, rational functions entries denoted with 

^ D + CiXIn-Ay'B, (4) 

Remark ILL Our results apply on both continuous or discrete time LTI systems, hence we 
assimilate the undeterminate A with the complex variables s or z appearing in the Laplace or 
Z-transform, respectively, depending on the type of the system. 

For elementary notions in linear systems theory, such as state equivalence, controlability, 
observability, detectability, we refer to [H], or any other standard text book on LTI systems. 

By M. pxm we denote the set ofpxm real matrices and by M(A) pxm we denote px m transfer 
function matrices (matrices having entries real-rational functions). 

This section contains a discussion based on reference ffT8l on the definition of the Dynamical 
Structure Functions associated with a LTI system. We start with the given system L(A) described 
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by the following state equations, of order n: 



x(t) = Ax(t) + Bu{t); x{t ) = x 
y(t) = Cx(t) 



(5a) 
(5b) 



Assumption IL2. (Regularity) We make the assumption that the C matrix from (TJQ) has full 
row rank (it is surjective). 



We choose any matrix C such that T 



def 



C 

c 



is nonsingular (note that such C always 



exists because C has full row rank) and apply a state-equivalence transformation 

x(t)=Tx(t), A = TAT~\ B = TB, C = CT' 1 . 
on (I5al).(l5bl) in order to get 
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Assumption II.3. (Observability) We can assume without any loss of generality that the pair 
(C, A) from d57il) . rt57il) or equivalently the pair (A12, A22) from ftTbl) are observable. 

Remark II.4. The argument that the observability assumption does not imply any loss of 
generality, is connected with the Leuenberger reduced order observer. 



Looking at the Laplace or Z-transform of the equation in (|7bl . we get 
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(note that 0(A) is always invertible as a TFM) we get 

((A/ p - An) - A 12 (A/ n _ p - Aaa)- 1 ^: 
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where the * denote entries whose exact expression is not needed now. Immediate calculations 
yield that the first block-row in ([TOl is equivalent with 



AF(A) = ( K A 11 + A 12 (XI n _ p ~A 22 )- 1 A 2i y(X)+^B 1 + A 12 (XI n _ p -A 22 )- 1 B 2 ) jU(X) (11) 
and by making the notation 



W(X) = f -A ll - A 12 (XI n - p - A 22 )- l A 2l 



V(X) d = f B 1 + A 12 (XI n ^-A 



p — ^22 j 



(12a) 
(12b) 



we finally get the following equation which describes the relationship between manifest variables 

XY(X) = W(X)Y(X) + V(X)U{X). (13) 

Remark U.5. Note that ifV(X) is identically zero, while ^(A) is a constant matrix having the 
sparisity of a graph's Laplacian, then ([73]) becomes the free evolution equation XY(X) = ^^(A). 
These types of equations have been extensively studied in cooperative control [15] to describe 
the dynamics of a large group of autonomous agents. Equation 4731) can be looked at as a 
generalization of that model and will be studied here in a different context. 

Since L(X) is the input-output operator from U(X) to Y(X), we can write equivalently that 
L(A) = (XI P -W{X))~ 1 V{X), which is exactly the (W, V) representation from Ml (3)/ pp.1671]. 
(Note that since W(X) is always proper it follows that (A/ p — W(X)) is always invertible as a 
TFM.) Next, let D(X) denote the TFM obtained by taking the diagonal entries of W(X), that 
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is D(X) = f diag{W u (A), W 22 (X) . . . W PP (X)}. Then we can write L(A) = (xi p - D(X) 
W(X) — D(X) ) y(X), or equivalently (note that (W — D) has zeros on the diagonal entries) 



L(A)= J- [XI P -D(X)) ( W(A) - -D(A) ) l (xi p -D(xj) 1 V(X) (14) 



and after introducing the notation 



Q(X) d = (xi p - D(X)^J (w{X)-D(X)) (15a) 
P(X)= f (aI p -£>(A))~V(A) (15b) 

we get that L(A) = (/ p - Q(X)j P(X) or equivalently that 

y(A) = Q(A)V(A) + P(A)C/(A) (16) 

Remark II.6. T/ze splitting and the "extraction " of the diagonal in A15a\) are made in order 
to make the Q(X) have the sparsity (and the meaning) of the adjacency matrix of the graph 
describing the causal relationships between the manifest variables Y(X). Consequently, Q(X) 
will always have zero entries on its diagonal. 

Definition II.7. / [7^1 Definition 1 ] Given the state-space realization ft7bl).ft7c\) of L(X) the Dynam- 
ical Structure Function of the system is defined to be the pair (Q(X), -P(A)), where Q(X), P(X) 
are given by U5a\) and A15b\) respectively. 



III. Dynamical Structure Functions Revisited 

One scope of this paper and also one of its contributions is to emphasize the idea that for 
a given TFM L(X) there exist more than one pair (Q(X), -P(A)) than the one in (|15al) . (115bl) 
(originally introduced in [18]) and which satisfy (fT6l . In fact there exists a whole class of pairs 
(Q(X), P(X)) that do satisfy (fT6l) and for which Q(X) has all its block-diagonal entries equal 
to zero. In order to illustrate this we need to slightly reformulate the original Definition 111.71 of 
Dynamical Structure Functions associated with a L(X) as follows: 
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Definition III.l. Given a TFM L(X), we define a Dynamical Structure Function representation 
of L(A) to be any two TFMs Q(X) E W pxp (X) and P(X) E R pxm (X) with Q(X) having zero 
entries on its diagonal, such that L(X) = (/ p — Q(X)) X -P(A) or equivalently 

Y{X) = Q(X)Y(X) + P(X)U{X) (17) 

The following definition will also be needed in the sequel. 

Definition III.2. Given the TFM L(X), we call a viable (W(X), V(X)) pair associated with L(X), 
any two TFMs W(X) E W xp (X) and V(X) E W xm (X), with W(X) having McMillan degree at 
most (n — p) and such that 



L(X) — (xi p — W(X)j V(A). (18) 

Proposition III.3. Given a TFM L(A) then for any given viable (W(X), V(X)) pair associated 
with L(X), there exists a unique DSF representation (Q(A),P(A)) of L(X) given by A15a\) and 
&15bk where D(X) d = diag{VFn(A), W^A) . . . W PP (X)} is uniquely determined by W(X). 

Proof: The proof follows immediately from the very definitions (I15al) . (I15bl) . ■ 

Remark III.4. It is important to remark here that any viable (^(A), V(X)) pair has the same 
sparsity pattern with its subsequent DSF representation (Q(A), -P(A)). For example W(X) is 
lower triangular if and only ifQ(X) is lower triangular. Similarly, for instance V(X) is tridiagonal 
if and only if P(X) is tridiagonal. 

Remark III.5. Using Proposition MIL 3\ we can conclude that in order to find all DSFs (according 
to Definition MIL 1 D associated with a given L(X), it is sufficient to study the set of all viable 
(w(X), V(X) j pairs associated with L(X). The following theorem gives closed-formulas for the 
parameterization of the class of all viable (w(X), V(X)^j pairs associated with a given TFM. 

Theorem III.6. Given a TFM L(X) having a state-space realization rf5a1).rt5/jl). we compute any 
equivalent realization rf7Z?l).rf7cl). The class of all viable (^W(X), V(X) j pairs associated with 
L(A) is then given by 
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where the K is any matrix in ]R(™~p) x p and An, A12, , A21, A22, B\, B 2 are as in d7^1).f[7cl). 

Proof: See Appendix B. ■ 

Remark III.7. We remark here the poles of both W(X) and V(X) can be allocated at will in 
the complex plane, by a suitable choice of the matrix K and the assumed observability of the 
pair (Ai2,A 2 2) (Assumption 1/7. 3D . 

IV. Main Results 

The ultimate goal of this line of research would be computing controllers whose DSF has a 
certain structure. This would allow us for instance to compute controllers that can be implemented 
as a "ring" network (see Figure [T) or as a "line" network which is important for motion 
control of vehicles moving in a platoon formation. However, classical results in LTI systems 
control theory, such as the celebrated Youla parameterization (or its equivalent formulations) 
render the expression of the stabilizing controller as a stable coprime factorization of its transfer 
function. As a first step towards employing Youla-like methods for the synthesis of controllers 
featuring structured DSF, we need to understand the connections between the stable left coprime 
factorizations (of a given stabilizing controller) and its DSF representation. We address this 
problem in this section. 

A. A Result on Coprimeness 

In this subsection we prove that (by chance rather than by design) for any viable (W(X), V(A)) 
pair associated with a given L(X) (with W(X) and V(X) as in Theorem IIII.6I) it follows that 
(Alp — W(X)) , V(X)) is a left coprime factorization of L(X). An equivalent condition for (XI p — 
W(A), V(A)) to be left coprime is for the compound transfer function matrix 



(xi p - W(\j) V(X) 



(21) 
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to have no (finite or infinite) Smith zeros (see |]3), DID, lITOl for equivalent characterizations of 
left coprimeness). Coprimeness is especially important for output feedback stabilization, since 
classical results such as the celebrated Youla parameterization, require a coprime factorization 
of the plant while also rendering coprime factors of the stabilizing controllers. 

Assumption IV. 1. (Controllability) From this point onward we assume that the realization 
rt5al).rf57il) of L(X) is controllable. 

Theorem IV.2. Given a TFM L{X), then for any viable (W(X), V(X)) pair associated with a 
given L{X) (with W(X) and V(X) as in Theorem \III.6\ it follows that (XI p — W(X), V^(A)) 15 a 
left coprime factorization of L(X). 

Proof: See Appendix B. ■ 

Remark IV.3. We remark here that while any viable (W{X),V{X)) pair associated with a 
given T(X) makes out for a left coprime factorization L(X) = (XI P — W(X)) 1 V(X), the DSF 
L(A) = (l p — Q(A)) 1 P(X) are in general never coprime (unless the plant is stable or diagonal). 
That is due to the fact that in general not all the unstable zeros of {XI p — D{X)j cancel out 
when forming the products in A15a\) . U5b\) and the same unstable zeros will result in poles/zeros 
cancelations when forming the product L(X) = {l p — Q(X)) X P{X). 

B. Getting from DSFs to Stable Left Coprime Factorizations 

In this subsection we show that for any viable pair (W(X), V(X)) with both W(X) and V(X), 
respectively being stable, there exists a class of stable left coprime factorizations. Furthermore, 
there exists a class of stable left coprime factorizations that preserve the sparsity pattern of the 
original viable pair {W(X), V'(A)). 

Note that for any viable pair {W(X), ^(A)) is an improper rational function and it has exactly 
p poles at infinity of multiplicity one, hence the {XI p — W{X)^ factor (the denominator of 
the factorization) is inherently unstable (in either continuous or discrete-time domains). We 
remind the reader that any the poles of both W(X) and V(X) can be allocated at will in 
the stability domain (Remark IIII.7I) . In this subsection, we show how to get from viable pair 
(W(X), V(X)) of L(A) in which both factors W(X) and V^(A) are stable, to a stable left coprime 
factorization L(X) = M -1 (A)iV(A). We achieve this without altering any of the stable poles of 
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W(X) and V(X) (which are the modes of (A 22 + KA 12 ) in ([191), d2Q])) and while at the same 
time keeping the McMillan degree to the minimum. The problem is to displace the p poles 
at infinity (of multiplicity one) from the (A/ p — W(X)) factor. To this end we will use the 
Basic Pole Displacement Result from |[T0l Theorem 3.1] that shows that this can be achieved by 
premultiplication with an adequately chosen invertible factor 6(A) such that when forming the 
product 6(A) (A/ p — W(X)) all the p poles at infinity of the factor (A/ p — ^(A)) cancel out. 
Here follows the precise statement: 

Lemma IV.4. Given a viable pair [XI p — W(X), V^(A)) of L(X) then for any 

©(A) 



def 


A x — XIp 


T 4 




T 5 






(22) 

with A X ,T^, T 5 arbitrarily chosen such that A x has only stable eigenvalues and both T4, T 5 are 
invertible, it follows that 



M(A) N(X) ] = 6(A) [ (xi p - W(\j) V{X) 
is a stable left coprime factorization L(X) = M~ 1 (X)N(X). Furthermore, 



(23) 
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(A x Ti - T A A U + T A A 12 K) T±B X 
(A22K + KA 12 K-KA 11 -A2i) KB 1 + B 2 




1 4 







I 



(24) 

hence all the modes in (A 22 + KAi 2 ) (which are the original stable poles ofW(X) and V(X)) 
are preserved in the M(X) and N(X) factors. 

Proof: See Appendix B. ■ 

Remark IV.5. We remark that for any diagonal A x having only stable eigenvalues 0(A) = 
{XIp— Az)^ 1 yields a stable left coprime factorization ofL(X) that preserves the sparsity structure 
of the initial viable ( AJ p — W(X), V(X) ) pair. 



C. Connections with the Nett & Jacobson Formulas KfflJ 

In this subsection, we are interested in connecting the expression from (|24l) for the pair 
(M(A), N(X)) to the classical result of state-space derivation of left coprime factorizations of a 
given plant originally presented in lfT6l (and generalized in [Tl7ll ). 
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Proposition IV.6. / I73I/ , / I77I/ Let L(X) be an arbitrary m x p TFM and Q a domain in C. The 
class of all /e/j coprime factorizations of L(X) over Q, T(A) = M _1 (A)iV(A), given Z?;y 



u- 1 



1 


(A - FC) - A/ 


-F 5 




C 


/ 



M(A) iV(A) 



where A, B,C,F and U are real matrices accordingly dimensioned such that 

i) U is any p x p invertible matrix, 

ii) F is any feedback matrix that allocates the observable modes of the (C, A) pair to Q, 



(25) 



Hi) L(X) 



A -XI 


B 


C 






is a stabilizable realization. 



Due to Assumption IIV.ll we have to replace the stabilizability from point Hi) with a contro- 
lability assumption. We start off with L(A) given by the equations (|5al) .(l5bl) 



L(A) 



A n - XI P 



A 



12 



A 21 A 22 - XI, 
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B 2 



O 



(26) 



and we want to retrieve (1241 ) by_ using the parameterization in Proposition IIV.61 First apply a 

T 4 O 



state-equivalence T 



K I 



in order to get 





n(A n - A 12 K)T A l - XI P 


T A A 12 


T 4 Bt 


L(A) = 


(KA U + A 21 - KA 12 K - A 22 K)T^ 


KA 12 + A 22 - XI n _ p 


KB X + B 2 




rp-l 









(27) 



Next, we only need to identify the F feedback matrix from point ii) of Proposition |IV.6[ which 
in this case is proven to be given by 



(T^A X T 4 ) - A u + A 12 K 
-K(T^A X T 4 ) + A 22 K - A 21 



(28) 



To check, simply plug (1281) in (|25l) for the realization (??) of L(X). 
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D. Getting from the Stable Left Coprime Factorization to the DSFs 

In this subsection we show that for almost every stable left coprime factorization of a given 
LTI system, there is an associated a unique viable (W(X),V(X)) pair and consequently (via 
Remark lni3l) a unique DSF representation (Q(X), P(X)). The key role in establishing this one 
to one correspondence is played by a non-symmetric Riccati equation, whose solution existence 
is a generic property. This result is meaningful, since for controller synthesis while we are 
interested in the DSF of the controller, in general we only have access to a stable left coprime 
of the controller. 

We start with a given stable left coprime factorization (|25) for L(X) having an order n 
realization 



M(A) N(X) 



1 


(A - FC) - XI 


-F B 




C 


I 



(29) 



to which we apply a type © state-equivalence transformation with T 6 IR nxn such that CT = 
I O ■ Note that such a T always exists because of Assumption 111.21 It follows that (1291 
takes the form 







A u + Ft- XI P 


A 12 


Fx 


Si 


M(A) N(X) 




A 2 i + F 2 


A 2 2 — AZ„_ P 


F 2 


B 2 

















and denote 



(30) 



A + 



(3D 



def ^-11 + -^1 —A12 

-(A 21 + F 2 ) A 22 

The solution to the following nonsymmetric algebraic Riccati matrix equation is paramount 
to the main result of this subsection, since it underlines the one to one correspondence between 
(1301) and its unique associated viable (W(X),V(X)) pair. 

Proposition IV.7. The nonsymmetric algebraic Riccati matrix equation 

K(A n + Fx) - A 22 K - KA 12 K + (A 21 + F 2 ) = O (32) 
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has a stabilizing solution K (i.e. (An + F x — A 12 K) is stable) if and only if the A + matrix from 
(1273 has a stable invariant subspace of dimension p with basis matrix 



Vl 

having V\ invertible (i.e. disconjugate). In this case K 



(33) 



V 1 1 V 2 and it is the unique solution 



Proof: It follows from D22J. ■ 

Remark IV.8. Since in our case A + is stable, all its invariant subspaces are actually stable 
( including the whole space ). Therefore, the Riccati equation has a stabilizing solution if and 
only if the matrix A + has an invariant subspace of dimension p which is disconjugate. Hence, 
if for example A + has only simple eigenvalues, the Riccati equation always has a solution (we 
can always select p eigenvectors (from the n eigenvectors) to form a disconjugate invariant 
subspace). In this case, all we have to do is to order the eigenvalues in a Schur form such 
that the corresponding invariant subspace has V\ invertible. Although this is a generic property, 
when having Jordan blocks of dimension greater than one it might happen that the matrix A + 
has no disconjugate invariant subspace of appropriate dimension p, and therefore the Riccati 
equation has no solution (stable or otherwise). 

Theorem IV.9. Given any stable left coprime factorization L(X) = M^ 1 (X)N(X) and its state- 
space realization (TJOj), let K be the solution of the nonsymmetric algebraic Riccati equation 021) 



and denote A 
given by 



def 



Fi + An — Ai 2 K). Then, a state-space realization for 



M(X) N(X) 



is 







A x AI n —p 


A u 


(A x - Au+A 12 K) Bi 


M (A) N(X) 




o 


A 22 + KA X2 - XI P 


(A 22 K + KA 12 K-KA 11 -A 21 ) KB 1 + B 2 






i 


O 


I 



(34) 



Furthermore, from ( TJ4l) we can recover the exact expression of the subsequent viable \ W( X) , V^(A) 
pair associated with L(X), where W(X) and V(X) are given by 4791) and rf20l) . respectively. 
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Proof: For the proof, simply plug 

A x - A u + A 12 K 
-KA X + A 22 K - A 21 



Fx 


def 


F 2 





(35) 



into the expression of (|30l ) in order to obtain ([34V The rest of the proof follows from Lemma llV.41 
by taking T 4 to be equal with the identity matrix I p . ■ 

Remark IV. 10. We remark here that in general there is no correlation between the sparsity 
pattern of the stable left coprime ([£2]) we start with and its associated viable (W(X),V(\)) 
pair produced in Theorem \IV.9\ That is to say that the converse of the observation made in 
Remark \IV.5\ is not valid. This poses additional problems for controller synthesis, since it might 
happen to encounter stable left coprime factorizations that have no particular sparsity pattern 
(are dense TFMs) while their associated viable (W(X), V(A)) pair are sparse. This is due to the 
fact that in general, the A x matrix in Theorem \IV.9\ can be a dense matrix. One way to circumvent 
this problem would be to use a carefully adapted version of Youla 's parameterization in which 
the stable left coprime factorization to be replaced with a DSF description where both with 
(W(A), V(X)) factors are stable. This is the topic of our future investigation. 

V. Conclusions 

In this paper we have presented an exhaustive discussion on the intrinsic connections between 
the DSFs associated with a given transfer function and its left coprime factorizations. We have 
showed that rather than dealing directly with the DSF representation it is more beneficial to work 
on the so-called viable (W(\),V(\)) pairs associated with a given system. This theoretical 
results ultimately aim at a method of designing LTI controllers that can be implemented over a 
network with a pre-specified topology. We currently have sufficient conditions for the existence 
of such controllers but we miss the necessary conditions. While in general these conditions 
might be very hard to find, we expect to find such conditions for plants featuring special DSF 
structures. 
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Appendix A 

Definition V.l. A TFM L(X) is called improper if for at least one of its entries (which are real- 
rational functions), it holds that the degree of the numerator is strictly larger than the degree 
of the denominator. 

Proposition V.2. f/[5]/, /[7|/j Any improper (even polynomial) p x m rational matrix L(X) with 
coefficients in IR has a descriptor realization of the form 

L(X) = D + C(XE- A)- 1 B =: 

where A,E G lR nxn , B G IR nxm , C G R pxn , D G R pxm , and the so called pole pencil A-XE 
is regular, i.e., it is square and det(A — XE) ^ 0. The dimension n of the square matrices A 
and E is called the order of the realization rf3ol) . 

Definition V.3. The descriptor realization |£6|) of L(X) is called minimal if its order is as small 
as possible among all realizations of this kind. 

Definition V.4. The McMillan degree of L(X) - denoted 8{L) - is the sum of the orders of all 
the poles of L(X) (finite and infinite). 

Remark V.5. The principal inconvenience of realizations of the form rfiol) is that their minimal 
possible order is greater than the McMillan degree of L(X), unless L(X) is proper, and this 
brings important technical difficulties in factorization problems in which the McMillan degree 
plays a paramount role. A remedy to this is to use a generalization of (TJH) in which either the 
"B" or the "C" matrix is replaced by a matrix pencil, as stated in the next Proposition. 

Proposition V.6. f/ fiOl/ ) Any improper p x m TFM L(X) has a realization 

= f D + C{XE - A)-\B - XF), (37) 
and for any fixed a, (3 G R, not both zero, there exists a realization 

W D + C(XE-A)- 1 B(a-Xf3), (38) 



A-XE 


B 


C 


D 



(36) 



L(A) 



A-XE 


B-XF 


C 


D 



L(X) 



A-XE 


B(a - A/3) 


C 


D 
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where A, E e R nxn , B,F G IR nxm , C G W xn , D G IR pxm , and the pole pencil A - XE 
is regular. A realization 081) will be called centered at % (if (3 = we interpret § as oo). 
Occasionally, we shall use also the more compact notation L(X) = (A — XE, B — XF, C, D) and 
L(X) — (A — XE,B(a — X(3),C,D) to denote 071) and 0ffl ), respectively. Realizations of type 
s ?!) have been dubbed pencil realizations. 

Definition V.7. f/fTOl/) We call realizations of the type 071) or s ?!) minimal if the dimension of 
the square matrices A and E ( also called the order of the realization ) is as small as possible 
among all realizations of the respective kind. 

Proposition V.8. Any TFM L(X) has a minimal realization of type 071 ) of order equal 

to 5(L). For any fixed a and (3, not both zero, and such that % is not a pole of L(X) there 
also exists a minimal realization of type 0ffl) of order equal to 5(L). The condition imposed on 
% is needed only for writing down minimal realizations 0ffl) which have order equal to S(L). 
More precisely, even if % is a pole of L(X) we can still write a realization 0^1) but the minimal 
order will with necessity be greater than 5(L). This is exactly what is happening for realizations 
0^1 ) which are obtained from 0^1 ) for a = 1 and (3 = 0, and for which the minimal order is 
necessary greater than S(L), provided § = oo is a pole of L(X). Notice that for 0g|) we can 
always choose freely a and (3 such as to ensure % is not a pole of L(X). For the rest of the 
paper, if not otherwise stated, we assume this choice implicitly. The nice feature of 071) and 0^1) 
that their minimal order equals the McMillan degree of L(X) recommends them for the kind of 
problems treated in this paper. 

Proposition V.9. A given realization of type 071) of a TFM L(X) is minimal if and only 

if all of the following conditions hold true 



rank A - XE B - XF 



n, VA G C 



(39a) 



rank E F 



n 



(39b) 



A — XE 



rank 



n, VA G C 



(39c) 



C 



E 



rank 



77 



(39d) 



C 
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while for realizations of type fi38\) similar conditions result by simply replacing (a) and (b) in 
(EB with 



rank 



A — XE B{a-\0) 



n, VA e C, 



rank 



E B 



n. 



(40a) 
(40b) 



Proposition V.10. Any two minimal realizations L(X) = (A — XE, B(a — X/3), C, D) and L(X) 
(A — XE, B(a — A/3), C, D) are always related by an equivalence transformation as 



E = QEZ, A = QAZ, B = QB, C = CZ, D = D, 
where Q and Z are unique invertible matrices. 



(41) 



Appendix B 

Proof of Theorem lHI.61 We prove that any pair (w(X),V(X)j given by (JT9]) , (J20l) satisfies 



181) . We start with the equations 



XI P - A u -A 12 




Y(X) 






—A 2 \ XI n ^ p — A 2 2 




. Z(A) . 






Y(X) =[l p o\ 




" Y(X) ' 
Z(A) 







U(X) 



and apply a type © state equivalence transformation with 



T 



I P O 



K L 



n—p 



where K can be any matrix in MS n p ) xp , in order to get 



(42a) 



(42b) 



(43) 



XI - (A n - A 12 K) -A 12 
(-KA n - A 21 + A 22 K + KA 12 K) XI - (A 22 + KA 12 ) 



Y(X) 



I P O 



Y(X) 
KY(X) + Z(X) 

Y(X) 
KY(X) + Z(X) 





Si 




K Bi + B 2 



U(X) 
(44a) 

(44b) 
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respectively. In a similar manner with getting from © to <TT0b via ®, we multiply (|44ab to the 
left with the following invertible factor 



i 



V 

o 



A 12 (XI, 



L n—p 



(45) 



- (A 22 + KA 12 ) 

In—p 

After the multiplication is performed, the first block row of the resulting equation yields (XI p — 
-W(X))Y(X) = V(X)U(X) which is exactly (JT8l) with W(X) and V(X) having the expressions 
in (TT9l) and (1201) . respectively. Finally, from the expression of W(X) in (fl9l) . clearly the McMillan 
degree of W(X) cannot exceed (n — p). 

Proof of Theorem IIV.2I An equivalent condition for the pair (A/ p — W(X),V(X)) to be 
coprime (over the compactification of C) is for the compound transfer function matrix 



XL 



(46) 



W(X)j V(X) 

to have no (finite or infinite) Smith zeros (see [0, flU, ifTOI for equivalent characterizations of 
left coprimeness). According to [10, Theorem 2.1] (see also 0, 0) the Smith zeros of (|46l are 
among the Smith zeros (generalized eigenvalues) of the system-pencil of any minimal realization 
of (l46l) . Hence we break this proof in two distinct parts: in part I) we compute a type (1381) pencil 
realization for (1461 ) and prove that is indeed minimal, in the sense of Definition IV. 9 1 In part II) 
of the proof we show that the system-pencil of the minimal realization from part I) has no finite 
of infinite Smith zeros (generalized eigenvalues). 

I) We will show that the following type (1381) pencil realization for [XI — W(X)) V(X) 
is a minimal realization in the sense of Definition IV.9I 





(A 2 2 + KA 12 )-XIn-p 


O 


{A 22 K + KA 12 K - 


KAu - An) KBi + B 2 


(XI-W(X)) V(X) 







I P 


Ip(Xo — 


A) 




A 12 


Ip 


X I P — An - 


f A i2 K Bi 



(47) 



I a) Observability for any finite X G C We note that 



(A 22 + KAu) - A/ n _ p 







I 





K 





h 







I 


-K 


A 12 












I 



A 22 - XL, 



n—p 



o 

A 12 



O 
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where the right hand side has full column rank for any A G C, due to the observability of 
the pair (Ai 2 ,A 22 ) (from Assumption 111.31 ) . Hence point (|39cb of Definition IV.9I holds via the 
Popov-Belevitch-Hautus (PBH) criterion. 

I b) Observability at X = oo is equivalent via point (|39dl) of Definition IV. 9 1 with the following 
matrix having full column rank 

/ O 
O O ■ 

A 12 I 

I c) Controllability for any finite A G C We look at the following succession of equivalent 
singular matrix pencils 



{A 22 + KA 12 ) - A/ n _ p O (A 22 K + KA 12 K-KA n -A 21 ) KB 1 + B 2 



O 



I P (Xo ~ A) 



O 



{A 22 + KA 12 )-XIn- P O {XK-KAn-An) KB l + B 2 



O 



UK - A) 



o 



{A 22 + KA l2 )-\I n - P K (X K-KA n -A 21 ) KB 1 + B 2 



O 



UXo - x) 



o 



(A 22 + KA 12 )-\I n _ p K (-KA 11 -A 21 ) KB 1 + B 2 



O 



-XL 



O 



{A 22 + KA 12 ) - A/„_ p K -A 21 KB l + B 2 
O I p -XI p + Au O 

{A 22 + KA l2 )-XI n . p K -A 21 B 2 
O I p -XI p + A n -B 1 



A 22 — A/ n _ p 
-A 12 



K 



I p —XIp + An 



-A 2 i 



B 2 
-Bi 



XI n - Ai 
-A 2X 



-A, 



XL 



n—p 



Bi —Ip 
A 22 B 2 K 



(48) 
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The full row rank of the last pencil above for any A G C, follows from the controlability 
Assumption IIV. 1 1 and the PBH criterion and it fulfills point (|40al) of Definition IV.91 
I d) Controlability at A = oo: is equivalent via point (I40bl) of Definition IV.9I with the following 
matrix having full row rank 

r j p (.;• o o 

O O I n - P o 

II) We look at the system-pencil of the realization (1471 ), namely 



S{\) d ^ 



(A 22 + KA 12 )-XI n - p O (A 22 K + KA 12 K-KA 11 -A 21 ) KB X + B 2 
O I p I P [Xo-X) O 



Al2 



X I p -A n + A 12 K 



(49) 



We will show next that the singular pencil in (1491 has no finite or infinite Smith zeros (generalized 
eigenvalues), which will conclude that the pair (xi — W{X),V{X) s j is left coprime. We will 
show this, by proving that S(X) keeps full row rank for any A G C and also for A = oo. 
II a) No Finite Smith Zeros We look at the following succession of equivalent matrix pencils 

(A 22 + KA 12 ) - XI n - p O (A 22 K + KA 12 K-KA 11 -A 21 ) KB 1 + B 2 



O I p 
A 12 I p 

{A 22 + KA 12 ) - A/„_ p 
O 

A 12 



I P (\o - A) 
A / p - A n + A 12 K 

O (\K — KA U — A 21 ) 
I p I p (X — A) 

I p X Ip — An 



O 
Bi 



A 



22 — A/ n _p 

O 

A 12 



Aoo - XI„ 



O -A 21 B 2 
I p I P (X -X) o 
Ip X I p — An Bi 



A 



\ 22 — aj„_p O —A 2 i B 2 
A l2 O XIp-Aii Bi 

O Ip Ip(X -X) o 



22 — A/„_p 
A l2 
O 



An - XI P 
A 2 i 
XL 



KBi + B 2 
O 
Bi 



O -A 21 B 2 
Ip X I p — An Bi 
I p I P (X -X) o 



A X2 
A 22 — A/ n _p 
O 



Bi O 
B 2 O 
O L 



The last pencil above clearly holds full row rank for any A G C due to Assumption IIV. 1 1 and the 
PBH criterion. 
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II b) No Smith Zeros at Infinity: Follows by the adaptation of Lemma 1]. 

Proof of Lemma IIV.4I This proof is based entirely on ffTOl Theorem 3.1] (Basic Pole Dis- 
placement Result). We start with the following type (|37l) minimal realization of 





Ip 





ip(A -A) 





(\I P -W(X)) V(X) 




O (A 22 - 


\- KA n ) - A/„_ p 


(A 22 K + KA l2 K - KA n - A 21 ) 


KBi + B 2 






Al2 


\ I P - An + A 12 K 


Bi 



(50) 

It can be observed that (|50l ) is already in the ordered block-Schur form [fTUl (2.14)/pp. 252]. We 
want to employ [10, Theorem 3.1] in order to compute the invertible TFM from [W[ (3.1)/pp. 
252] which we denote with 0(A) that by premultiplying (|50l will cancel out the p poles at 
infinity of (|50l ). Any type (1381) realization of a valid 6(A) satisfies [fTOl (3.2)/pp. 252] for certain 
invertible X and Y matrices. Hence for any 



©(A) = 



A x - XI P 


B X (X - X ) 


C x 





(with D x must be invertible because 6(A) is invertible ) we write the conditions from [[TO] 
(3.2)/pp. 252] which are equivalent with 

A x — XIp B X (X — X Q ) 
C x D x 

From the first block row of (1511 we get that C X X = —D x and from the second block-row of 
(I5TT) we get B X (X - X Q ) = Y - (A x - XI P )X. Consequently 



©(A) = 



A x — XIp 


Y - (A x - XI P )X 


C x 


—C X X 



which is equivalent with 



©(A) = 



A x — XIp 


Y 


c x 






where C x and Y are arbitrary invertible matrices. We have denoted C x with T 4 and we have 
denote Y with T 5 to avoid notational confusion. The proof ends. 





X 




Y 




I 








{I P - AO) 



(51) 
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